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Abstract 

Different features of a potential in the form of a Gaussian well have been discussed extensively. 
Although the details of the calculation are involved, the general approach uses a variational method 
and WKB approximation, techniques which should be familiar to advanced undergraduates. A 
numerical solution of the Schrodinger equation through diagonalization has been developed in a 
self-contained way, and physical applications of the potential are mentioned. 



1 



I. INTRODUCTION 



Potentials such as the infinite square well, the harmonic oscillator, the delta function 
well, and the finite square well are frequently discussed in textbooks^ as examples that 
have bound states. In this paper we consider solutions of the time independent Schrodinger 
equation^ 

- 2^^0*0 + V(x)^(x) = (1) 
for the Gaussian well given by 

V(x) = -V e~ ax2 (-00 < x < +00), (2) 

where Vq > and a > 0. In Eq. (CEJ) m is the mass of the particle and ip(x) is the eigen- 
function. We will obtain an estimate of the ground state energy from a simple variational 
method and determine that there are a finite number of bound states using the WKB ap- 
proximation. We also formulate Eq. ([1]) as a matrix eigenvalue problem, which can be used 
for any Hamiltonian. For simplicity, we choose units such that h — m — 1 in the following. 



II. BOUND STATE CRITERION 



Consider the Hamiltonian H for a particle in one dimension subjected to a potential 

V(x), 

*=-i£ + V(,). (3) 

To demonstrate that H possesses a bound state, it is sufficient to construct a real normalized 
trial function ip(x) such that the expectation value of H for this trial function is negative, 
that is, 

r+00 

(H) = / ipHipdx < 0. (4) 



This sufficiency condition becomes 

{H) = 



dx. 



(5) 



Defining L to be the width (or, more precisely the length scale) of ip(x), we can write ip in 
the form 



We substitute this form for ip in Eq. (jSJ) and obtain 

(H) 



ai 02 

2T2 + T : 



(7) 



where a\ = J_^(d(j)/dx') 2 dx' (with x' = x/L), is a dimensionless positive number and 
a<i = j_ 00 V(p 2 dx. If L is sufficiently large the second term on the right-hand side of Eq. ([7]) 
dominates the first term, and [0(x/L)] 2 becomes ~ [0(O)] 2 (a positive constant), so that, it 
can be taken out of the integral. Hence, the sign of a 2 depends only on the sign of Vdx. 
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For both a± and a-i positive, it is not possible to obtain a positive value of L for which (H) 
is negative. Thus, the condition on V(x) to have at least one bound state is 

/+00 
Vdx < 0, (8) 
-00 

which depends on the shape of the potential, not on its strength. 

Equation (jSj) is a sufficient condition for a potential to have a bound state. It is not a 
necessary condition as the example of the harmonic oscillator potential shows. Note that 
this general criterion for the existence of the bound state follows from a simple argument 
involving dimensional analysis and the variational principle. 



III. GROUND STATE ENERGY FROM THE VARIATIONAL METHOD 

In practice, we can obtain various upper bounds of the ground state energy by calculating 
(H) for suitably chosen trial functions. There is a tradeoff between improving the approx- 
imate ground state energy using a complex trial wave function and the ease of calculation. 
Here we choose a normalized Gaussian trial wave function ip(x) with adjustable width, 

fw-dyv* (9) 

where b is related to the width L. The expectation value of the Hamiltonian H becomes, 

(H) = --V 0X I-^. (10) 
v 1 2 V 2b + a K J 

From Eq. ( flOl) the condition for a minimum is d(H)/db = 0, or 

b(2b + a) 3 = 2V r V. (11) 



Vb 


a 


b 


(H) 


E 


1.0 


1.0 


0.3742 


-0.4671 


-0.4774 


2.5 


0.5 


0.6113 


-1.8005 


-1.8038 


3.0 


1.0 


0.8717 


-1.9557 


-1.9637 


3.0 


0.1 


0.3504 


-2.6312 


-2.6316 



TABLE I. Comparison of the ground state energy (H) calculated from the variational principle 
and the numerical solution (E num ). (H) is within 2% of the corresponding numerically evaluated 
value. 

Equation fllljl determines the variational parameter b for given values of Vq and a. The 
variational method ground state energies for different values of Vo and a are given in Table [B 
As expected, these estimates are greater than the exact (numerically obtained) ground state 
energies, but are remarkably close. 
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IV. FINITE NUMBER OF BOUND STATES 



The wave function ip n (x) corresponding to the eigenvalue E n of a discrete spectrum has 
n nodes. The nature of the number of bound states is best demonstrated by the WKB 
approximation. States belonging to the discrete energy spectrum are semi-classical only for 
large values of n. In one dimension the WKB integral for the energy E is given by2 



rx2 

/ y/2[E-V(x)]dx 



n 



- |77. 

2/ 



for n = 1, 2, 3 . . . 



(12) 



where X\ and x-i are the turning points of the classical motion. Because the discrete spectrum 
lies in the range of energy values for which the particle cannot move to infinity, the energy 
must be less than the limiting values V(x = ±oo). For the Gaussian potential this condition 
implies that E < 0. Thus, the number of discrete levels iV is obtained by setting E = in 
Eq. f lT2"j) . so that n = N is close to the quantum number for the last bound state. For the 
classical turning points x = ±oo, Eq. f[T2"j) becomes 



V 7 ^ 



o / e 

oo 



- ax2 ' 2 dx 



(13) 



which gives 

y/n V a 2 

For Vo and a finite, N is also finite. Therefore, the number of bound states is finite for 



N 



(14) 



Vo 
a 


JV 


^num 


0.5 


1.3 


1 


1.0 


1.6 


1 


10.0 


4.1 


4 


100.0 


11.8 


11 



TABLE II. Number of bound state given by the WKB approximation (iV) and the numerical 
method (A^" num ). 



the Gaussian well. The value of iV calculated from Eq. ( IT4l) compared to the numerical 
solution is in Table HI1 Equation ( IT4"1) implies that the number of bound states is a function 
of the ratio Vo/a, which is a measure of the length scale L of the wavef unctions.- Because 
V(x) oc (y^a/7r)e~ ax , V(x) becomes a delta function in the limit a — > oo, for which Eq. ( |T3|) 
gives iV = 1/2, implying that there exists exactly one bound state for the delta function 
well.- The reason for the half-quantum number is that the ground state has no nodes. If 
Vo — > oo, then iV — y oo, which is the case for an infinite square well.— 



V. NUMERICAL SOLUTION OF THE SCHRODINGER EQUATION 

The one-dimensional Schrodinger equation has few analytic solutions, and most problems 
must be solved numerically. We first use a Taylor's series to obtain a discretization of 
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derivatives for a function f(x) and write the second derivative of the function f(x) in the 
three point central difference form 



/"(*) 



f(x + 5)-2f(x) + f(x-6) 
8 2 



+ 0(5 2 



where, 5 is the step size. If r is the number of mesh points and x 
and minimum value of the variable x, the step size 5 is 

5 



max; -^min 



(15) 

are the maximum 
(16) 



The solution of Eq. (DO) is accomplished by discretizing it using Eq. ( f!5|) and evaluating 
the functions and derivatives at x k = x min + kd (for k — 1, 2, . . . , r — 1), 



j)(x k + 5) - 2j)(x k ) + jj(x k - 5) 
25 2 



+ V(x k )i)(x k ) = Ei[)(x k ), 



(17) 



Because the solution is expected to decay exponentially outside x max and a; mm , we solve 

]. These cut-offs and the step size need to be adjusted to 



Eq. (fTT|) in the interval [x max , x minJ 
obtain the desired accuracy. 

Equation (TTTj) is equivalent to the tight-binding approximation applied to a chain of atoms 
with spacing 5 and one orbital per atom.- We can develop a matrix representation of the 
Schrodinger equation, with —5~~ 2 /2 in the sub-diagonal and super- diagonal matrix elements 
and the diagonal elements are 5~ 2 + V(x k ). Equation (TTTj) can then be written as a matrix 
equation as, 



(b~ 2 + V( Xl ] 
-5- 2 /2 






-5- 2 /2 
S- 2 + V{x 2 ) 





-5- 2 /2 













ip(x 2 ) 



5- 2 + V(x r _ 2 ) -5- 2 /2 
-5- 2 /2 5- 2 + V(x r ^)J 



1p(x r - 2 

VV(a?r-i)/ 



E 



( I'M \ 

ip(x 2 ) 

ip(x r _ 2 

\lp(x r -!)J 



(18) 



Equation ( fl8l) is a matrix eigenvalue problem with a tridiagonal matrix of dimension 
(r — 1) x (r — 1); thus there are (r — 1) eigenvalues. Because all the matrix elements are 
real and the transpose of the matrix is equal to the matrix itself, it is a Hermitian matrix, 
and hence all the eigenvalues are real. An efficient way of diagonalizing a tridiagonal matrix 
is to use the standard LAPACK routine DSTEVX, 5 which stores the symmetric tridiagonal 
matrix in two one-dimensional arrays, one of length (r — 1) containing the diagonal elements, 
and one of length (r — 2) containing the off-diagonal elements, and returns the eigenvalues 
along with the eigenfunctions of the matrix. 

The eigenvalues and corresponding eigenfunctions for the Gaussian well with Vq = 3.0 
and a = 0.1 are shown in Fig. [lj The ground state is symmetric with respect to the center 
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FIG. 1. The energy eigenfunctions and the corresponding eigenvalues for V$ = 3.0 and a = 0.1. 
Note the unequal spacing between different levels. 

of the potential (even parity), and the first excited state is antisymmetric (odd parity). The 
wave functions resemble those of the harmonic oscillator,- because for x < 1, the dominant 
term in the expansion of e~ ax is proportional to x 2 . 

The virtue of the matrix method is that it can be applied to any potential for which the 
Hamiltonian can be brought into a symmetric tridiagonal or bidiagonal matrix. However, 
for systems with periodic boundary conditions, the Hamiltonian is no longer tridiagonal and 
we cannot use the simple matrix method. 



VI. TUNNELING IN A GAUSSIAN BARRIER 

A Gaussian barrier can be constructed by changing the sign of Vq in Eq. 02]). Because the 
condition in Eq. (jSJ) is not satisfied, there is no bound state. The barrier formed in this way 
is an interesting example of tunneling. We start with the time energy uncertainty principle 

AEAt ~ -. (19) 

Denote the energy of the incident particle by E. The uncertainty in the energy is AE and 
for sufficiently small At the energy of the particle E + AE is greater than the height of the 
barrier Vq. Tunneling will take place if in the time At the particle can traverse the barrier. 
We take a -1 / 2 as the width of the barrier and write, 

a" 1/2 , x 

At ~ . 20 

y/2(E + AE- V ) 
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From Eqs. ( !T9|) and ( f20]) we find AE to satisfy the equation, 

(A£) 2 -|A£ + |(Vo-£) 
The condition for A.E to be real implies that 



0. 



a 



>(Vo-E), 



(21) 



(22) 



which is the condition for tunneling. The left-hand side of Eq. (1221) is the kinetic energy of 
the particle obtained from the 'position-momentum uncertainty relation, with an uncertainty 
of Ax ~ a~ 1//2 in the particle's position. We see that for tunneling to occur, the kinetic 
energy of the particle must be greater than the difference between the height of the barrier 
Vo and the total energy E. 

We now investigate tunneling using the WKB approximation. Consider the classically 
inaccessible region, E < V(x), as a < x < b, so that we can write the transmission coefficient 



as 2 



T 



cxp 



n(x)dx 



where k(x) 



2(V(x) — E). We define the opacity of the barrier by 

r r b i 

9 = exp / n(x)dx , 



(23) 



(24) 



so that Eq. (|23|) becomes T « Q- 2 . No attempt been made to obtain an analytical solution 
to the integral in Eq. ([23]) . But we can find an approximate solution of the integral that 
correctly predicts the transmission coefficient T. 

We introduce the dimensionless quantity — Vq/E > 1 and write the classical turning 
points as a - 
Eq. ([21]) to 



^/\n/3/a and b = ^/\n/3/a. The change in variables, x = y^lnfi/a, changes 

01 _(!_£-»» + 0-1)]^ 
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cxp 



'2Vbln/3 



(25) 



The term (1 — y + (3 : ) is always less than unity. If we assume that this term is 



approximately equal to unity, 6 then the binomial expansion of 
first order we have 



9 ~ exp 

which can be further expressed as 
9 ~ exp 



2V \n/3 



1 



i - i - i3~y 2 + (3- 1 



dy 



to 



(26) 




VK5 + ^erf(VK5) - 



(27) 



where erf(x) is the error function of x. Equation ( [271) implies that T depends only on the 
ratio — Vq/E, a common feature of barrier tunneling. The dependence of T on /3 can be 
seen from Fig. |2j It is evident that the transmission coefficient decreases as f3 increases. As 
expected for — 1, that is, E = V(x = 0), T = 1, implying no reflection from the barrier. 
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FIG. 2. Dependence of the tansmission coefficient T on j3 = Vq/E. As expected, T decreases as /3 
increases. 

VII. TWO SIMPLE APPLICATION 

Thus far we have discussed different aspects of the potential from a quantum theoretic 
approach. The Gauusian well (or, barrier) is not a long range potential, because it falls off 
faster than \. However, it has a crucial advantage over the widely used 'finite square well' 
potential in that it is continuous throughout the entire range [— oo < x < +oo], whereas the 
finite square well goes to zero discontinuously. In the following we describe two applications 
of this potential. 

A. Single Particle Motion in Atomic Nuclei 

The 'mean-field' dynamics of a single nucleon in the field of all other nucleons is the 
starting point of nuclear many-body theory. The Gaussian well potential is well suited for 
describing the interaction of a nucleon (especially, a neutron) with the heavy nucleus.- The 
force between them is appreciable only over a very short distance, of the order of 10 _15 m, and 
in this range the forces are very large compared to forces holding atoms together. Keeping 
in mind the three dimensionality of a nucleus, we can change our variable from x to the 
radius vector r in spherical polar coordinates, so that the potential can be represented by a 
'half Gaussian well, 



It is monotonically increasing with distance (i.e. attracting) and approaches zero very quickly 
as r goes to infinity, reflecting the short-range nature of the nuclear force. The parameter 




r > 



r < 0. 



(28) 
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a can be adjusted to the experimental value of the nuclear interaction barrier. From the 
Schrodinger equation in radial form, one has 



2m dr 2 



+ 



■ v e 



2m 



u 



Eu, 



(29) 



where / is the azimuthal quantum number. It follows immediately that this equation is 
identical to Eq. p]) for I = 0. It is left as an exercise to the reader to solve this equation 
using the numerical method discussed already and to deduce the single particle nuclear 
energy levels for / = by taking appropriate values of Vq and a.— For a potential such as 
the one in Eq. ( )28l) there will be no even-parity states (because at r = there must be an 
antinode.). 



B. Scanning Tunneling Microscope 

The Gaussian barrier can also be very useful in analyzing a scanning tunneling microscope 
(STM). The basic idea behind a STM is quite simple. An atomically sharp metal tip is 
brought very close (< 5 A) but without physical contact to a sample surface, and a small 
bias voltage (3-6 volts) is applied between them. If the distance is small enough, an electron 
(with an energy 3-6 eV) can tunnel quantum mechanically through the potential barrier 
developed in between the tip and sample. This gives rise to the tunneling current that is the 
result of the overlap between wave functions of the tip-atom and the surface-atom. Assume 
the metal to tip gap is a Gaussian barrier with a height of Vq = 10 eV ~ 0.36 in atomic 
units. An electron with energy of E = 5 eV ~ 0.18 in atomic units approaches the surface. 
Now from Eq. ( 123]) . the transmission coefficient is 

T ~ e" 2 - 2 ^ (30) 

If the tip is = 0.15 nm ~ 2.83 in atomic units from the surface, we have T ~ 0.024. 
However, it should be noted that for a large distance between the tip and surface, e.g., if 
-5= = 1 nm ~ 18.9 in atomic units, T will become ~ 10 -11 and practically, there will be 
no current to measure ! Thus the magnitude of the tunneling current is extremely sensitive 
to the gap between the tip and sample. The position of the tip in three dimensions is 
accurately controlled by piezoelectric drivers. The tip is scanned in two lateral dimensions, 
while a feedback circuit constantly adjusts the tip height to keep the current constant. As 
we measure the current with the tip moving across the surface, information about the atomic 
nature of the surface can be determined by tracing the path of the tip. 

As a final example there is the alpha particle tunneling problem which also can be modeled 
with the Gaussian barrier. 



VIII. A DOUBLE GAUSSIAN WELL 

An interesting double well potential may be formed by multiplying V(x) in Eq. ([2]) by 
a factor of x 2 . Because the potential U(x) = —Vox 2 e~ ax is symmetric, the Hamiltonian 
can be brought into a symmetric tridiagonal form. The numerical solution for this potential 
is given in Fig. |3j Because the barrier is finite the wave function extends into both wells. 
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The ground state wave function is small but nonzero inside the barrier. Assume the particle 
starts in the right well. From Fig. [3] we see that initially ip — ipi (where, ipi & ^2 are the 
ground state and first excited state, respectively) is canceled on the left. As a consequence 
the particle will oscillate between the wells with the period, r = 27ih/(E 2 — E\). Thus the 
tunneling rate depends on the energy difference, AE = E 2 — E\, and a double well with a 
high or wide barrier will have a smaller AE than one with a low or narrow barrier. Also, 
AE will become larger as the energy increases (that is, as {V — E) decreases). 




-4-3-2-101234 

Position (a.u.) 



FIG. 3. The energy levels (ground state & first excited state) along with the wavefunctions for 
U(x) with Vq = 3.0 and a = 1.0. We can see the tunneling of the wavefunction through the barrier. 



The double Gaussian well is a good model for a two level system. An example is a 
quantum well laser, based on the transition that an electron makes between the ground 
and first excited state of a double Gaussian well. By choosing a and Vq, we can tune 
the wavelength of the light emitted. Another not so obvious case could be a Qubit. The 
quantum bit, or qubit, is the simplest unit of quantum information. Measurements give only 
two values: 'zero' or 'one'. It is described conveniently by a state in a two level quantum- 
mechanical system. A pure qubit state would then be a linear superposition of those two 
states. Though the double Gaussian well has multiple energy levels, it can be seen from 
Fig. flU) that the relative spacing between the ground and first excited state is very small 
compared to that between the first and second excited states. Therefore, the lowest two 
bound states can be effectively decoupled from the other states by choosing appropriate 
values of Vq & a. An electron making transitions between these two states can easily be 
considered as a 'charge qubit', a much discussed way for consistent quantum data storage. 
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FIG. 4. The wavefunctions and the corresponding energy eigenvalues of U(x) for Vq = 10.0 and 
a = 1.0. The proximity between the lowest two states is evident. 

ACKNOWLEDGMENTS 

The author would like to thank Sayan Chakraborti, Swastik Bhattacharya, and Shamashis 
Sengupta for many helpful suggestions. He is grateful to Professor Deepak Dhar for help in 
understanding different aspects of quantum mechanics. 



* saikat@tifr.res.in 

1 D. J. Griffiths, Introduction to Quantum Mechanics, 2nd ed. (Prentice Hall, NJ, 2005), Chap. 2. 

2 E. Merzbacher, Quantum Mechanics, 3rd ed. (John Wiley & Sons, NY, 2005), Chap. 7. 

3 Here v/l/a is the length scale for the potential. Now keeping the ratio (Vo/a) fixed, the length 
scale of the wavefunctions can be identified as L ~ ~y^j^ = v^/^o- 

4 T. B. Boykin, "Tight-binding-like expressions for the continuous-space electromagnetic coupling 
Hamiltonian," Am. J. Phys. 69, 793-798 (2001). 

5 E. Anderson et al., LAPACK Users' Guide, 3rd ed. (SIAM, Philadelphia, 1999), p. 338. 

6 We can add more terms in the binomial expansion to obtain more accurate values of T. However, 
the essential physical behavior of T will remain the same. 

7 It can be thought as a potential similar to the well known Woods-Saxon potential. The 
Schrodinger equation for Woods-Saxon potential has been solved for / = to obtain the nu- 
clear energy levels. See for example, S. Fliigge, Practical Quantum Mechanics, 2nd ed. (Springer 
Verlag, Berlin, 1999), pp. 162-166. 

8 In principle, one can also solve for / ^ by this numerical method. The diagonal terms in the 
tridiagonal matrix will then be modified accordingly. 



11 



